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Abstract 

According to first-principles calculations, the solubility of Mg as a substitute for Ga or Al in Al x Ga-| _ X N bulk is limited 
by large, positive formation enthalpies. In contrast to the bulk case, the formation enthalpies become negative on 
Al x Ga-| _ X N surface. In addition, the N-rich growth atmosphere can also be favorable to Mg incorporation on the 
surface by changing the chemical potentials. On the basis of these special features, we proposed a modified surface 
engineering technique that applies periodical interruptions under an ultimate V/lll ratio condition (extremely N-rich), 
to enhance Mg effective incorporation. By optimizing the interruption conditions (2 nm interruption interval with 
2 s interruption time), the enhancement ratio can be up to about 5 in the AI 099 Gao.oiN epilayer. 
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Background 

Al^Gai _ alloys have attracted considerable attention 
in recent years because of their great potential for applica- 
tions in UV and deep UV optoelectronic devices with 
spectral lengths as short as 200 nm [1]. Both high-quality 
p-type and n-type AIGaN epilayers are strongly demanded 
for electrical injection in constructing these short wave- 
length devices. However, similar to most wide bandgap 
semiconductors, AIGaN suffers from the asymmetric dop- 
ing' limitation [2,3], i.e., doping AIGaN to form w-type 
layer is easy, but achieving p-type doping is difficult [4,5]. 
Although Mg is the most widely adopted p-type dopant 
for AIGaN, its doping efficiency is extremely low, particu- 
larly for high Al content Al^Ga! _ X N [6]. The low doping 
efficiency of Mg is mainly attributed to its limited solubil- 
ity, high activation energy, and compensation effect with 
impurities or native donor defects [2,7]. In spite of the 
extensive efforts to improve the Mg activation efficiency 
[5,6,8,9], the bottleneck of low Mg solubility in GaN [10] 
and A1N [11] materials strongly restricts the overall p-type 
doping in AIGaN. 
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Regarding the dopant solubility issue, an extremely 
high carbon dopant concentration was shown to exist 
on the epitaxial surface of Si system [12], This high con- 
centration can be attributed to the surface enhancement 
effect caused by the partial release of atom mismatch 
strain. As the epitaxy continues, part of this high con- 
centration dopant segregates to the new surface, and the 
residual components freezes into the host matrix [12] 
which corresponds to the final dopant concentration. In 
other words, the growing surface plays a critical role in 
determining dopant solubility. If the transient solubility 
on the growing front surface can be effectively enhanced, 
high dopant incorporation can be achieved. Theoretical 
simulations have recently predicted that a N-rich con- 
dition is beneficial for Mg incorporation in GaN and 
A1N [10,11]. However, high V/III ratio was determined 
to be unfavorable for high-quality Al^Ga! _ ^N crystal 
growth [13-16]. Thus, the dilemma between maintain- 
ing high V/III ratio to promote Mg incorporation and 
maintaining low V/III ratio to ensure high crystal quality 
presents a long-standing challenge for deep UV optoelec- 
tronic devices. 

In this work, we proposed a method to solve this V/III 
ratio dilemma by periodically interrupting the AIGaN 
growth (using usual V/III ratio as the AIGaN growth) 
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and by shortly producing an ultimate V/III ratio condition 
(extremely N-rich). First-principles simulations were 
utilized to analyze the behavior of substituting Mg for 
Al and Ga in the bulk and on the surface of Al^Ga! _ 
under different growth atmospheres and to demonstrate 
the mechanism for the preferred Mg incorporation. On 
the basis of the analysis results, a modified surface engin- 
eering (MSE) technique that utilizes periodical interrup- 
tions under an extremely N-rich atmosphere was applied 
to enhance Mg effective incorporation by metalorganic 
vapor phase epitaxy (MOVPE). Significant Mg incorpor- 
ation improvements in Al-rich Al^Ga! _ JST epilayer were 
achieved. 



Methods 

The first-principles total energy calculations based on 
density functional theory were performed by using the 
Vienna ab initio simulation package [17]. Pseudopoten- 
tials were specified by the projector augmented wave 
[18,19] and by generalized gradient approximation [20]. 
Ga 3d electrons were treated as part of the valence 
band, and the plane wave cutoff energy was set at 
520 eV. Geometry optimizations were performed until 
the total energy converged to 1 meV. For the bulk calcu- 
lations, a 2 x 2 x 4 supercell containing 64 atoms [7] and 
a 5x5x3 Monkhorst-Pack grid [21] of /c-points were 
used. All atoms were allowed to relax fully for energy 
minimization. For the surface calculations, we employed 
a 2 x 2 supercell with six Al^Ga! _ bilayers separated 
by a 13-A wide vacuum region [22] and a 4x4x1 
/c-point mesh. The back side of the slab was saturated 
with hydrogen atoms of fractional charge. The three bot- 
tom Al^Ga! _ JST bilayers were fixed in the appropriate 
bulk-optimized configuration to simulate the growth sur- 
face, in which all the other layers was relaxed fully. 

The Mg-doped Al^Ga! _ samples were grown on 
(0001) sapphire substrates via MOVPE. Trimethylgallium 
(TMGa), trimethylaluminum (TMA1), bis-cyclopentadienyl- 
magnesium (Cp 2 Mg), and ammonia (NH 3 ) were used as 
precursors, and H 2 was used as carrier gas. Buffer layers 
with a 20-nm low temperature A1N nucleation layer, a 
1-um high temperature A1N layer, and a graded com- 
position AlGaN layer have been used for initial growth 
on sapphire. The conventional method for fabricating 
Mg-doped Al^Gai _ was conducted by retaining all 
the precursors during Mg-doped growth with usual V/III 
ratio. The MSE technique was implemented by period- 
ically interrupting the conventional growth mode with 
closing the metal flows (TMA1, TMGa, and Cp 2 Mg) 
and continuously maintaining the NH 3 flow to shortly 
produce an ultimate V/III ratio. The Mg and H concentra- 
tions were measured by using the Quad PHI 6600 second- 
ary ion mass spectrometer (SIMS) system with depth 



resolution of approximately 2 nm, and Cs + ion beams were 
used as primary ion sources. 

Results and discussion 

Considering that MOVPE growth is usually characterized 
by N-rich growth, we first discuss the formation enthalpies 
of neutral charge state Mg substituting for Al (Mg A i) and 
Ga (Mg Ga ) in Al^Ga! _ *N bulk as a function of Al content 
under N-rich condition. The calculated results are shown 
in Figure la, wherein both the Mg A i and Mg Ga formation 
enthalpies are positive and large, thus indicating limited 
Mg solubility. The formation enthalpies of Mg A i in A1N 
and Mg Ga i n GaN are comparable with previous results 
[10,11]. As the Al content in Al^Ga! _ ^N increases, both 
the Mg A i and Mg Ga formation enthalpies monotonically 
increase. The formation enthalpy AHf is closely related to 
the equilibrium Mg solubility C, which is given by [10]: 

C = N sites e~ AH f^ T (1) 

where N sites is the number of sites on which the dopant 
can be incorporated, k B is the Boltzmann constant, and 
T denotes the temperature. Large formation enthalpy 
yields low dopant solubility. At the growth temperature 
(T= 1,000°C), the Mg solubility in bulk GaN is approxi- 
mately 1.65 x 10 17 cnT 3 . Considering that AHf increases 
with increasing Al content, Al^Gai _ ^N experiences an ag- 
gravating Mg solubility limit. The Mg solubility limit may 
even decrease to approximately 2.32 x 10 16 cm" 3 in A1N 
(for T= 1,200°C). On the basis of this tendency, incorpor- 
ating Mg becomes more difficult in Al-rich Al^Ga! _ X N. 
Notably, the formation enthalpy for Mg A i is larger than 
that for Mg Ga over the entire Al content range. This char- 
acteristic demonstrates that substituting Mg for Al is more 
energetically unfavorable than substituting Mg for Ga, 
which also explains the low Mg incorporation in Al-rich 
Al^Ga! _ ^N. Such behavior of Mg is partly attributable to 
its larger covalent radius (1.36 A) compared with those of 
Al (1.18 A) and Ga (1.26 A), as well as the compressive 
strain after Mg substitution [23,24]. As shown in the 
inset of Figure la, the Al^Ga! _ lattice constants a 
and c decrease as the Al content increases, thus making 
the mismatch strain caused by substituting Mg for Al or 
Ga atoms with smaller radii becomes more considerable. 

Reducing the formation enthalpy is believed to be the 
key issue in solving the problem of Mg incorporation. 
The formation enthalpy is governed by two important 
factors, as given by [11]: 

AHf {neutral state) = AE + Aft. (2) 

Here, AE = E Mg - E host , where £ Mg and £ host are the 
total energies of the supercell with and without Mg sub- 
stitution; A(A = //Ai/Ga - {*Mg> where ^ (i = Al, Ga, Mg) is 
the chemical potential. AE is mainly induced by the 
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Figure 1 Formation enthalpies of Mg Ga /Mg A | and normalized 
Q/ig cprofile of AIGaN films, (a) In the bulk and (b) on the surface 
of Al x Ga] _ X N as a function of Al content under N-rich condition, 
(c) Normalized C Mg of Al x Ga-| _ X N {x = 033, 0.54) epilayers from the 
surface to bulk. The inset in (a) shows the calculated Al x Ga] _ X N 
lattice constants a and c as a function of Al content. The inset in 
(c) illustrates the source supply sequence of the conventional 
method. 



strain caused by the atom size mismatch. Consequently, 
larger atom size mismatch results in large A£, thus 
resulting in larger AHf as mentioned above. The strain 
induced by the C-dopant into the Si matrix becomes 
smaller on the surface than in the bulk [12]. The question 
of whether the surface also plays a similar role in the Mg 
incorporation in Al^Ga! _ arises. To address this issue, 
we further investigated the formation enthalpies of Mg Ai 



and Mg Ga on Al^Ga! _ surface, and the results are 
shown in Figure lb. In contrast to the bulk case, both of 
the formation enthalpies in the surface are negative, sug- 
gesting favorable Mg substitution. The value of Mg A i be- 
comes lower than that of Mg Ga and decreases as the Al 
content increases. These interesting reversed tendencies 
provide us a possible way to promote the Mg incorpor- 
ation in Al-rich Al^Ga! _ by utilizing the surface 
effect. An epitaxy growth, e.g., MOVPE and molecular 
beam epitaxy systems, is conducted under an inherently 
non-equilibrium process with the surface transforming 
into a bulk [12]. Therefore, enhancing the Mg incorpor- 
ation by using the surface effect should be practically 
feasible. 

Two Mg-doped Al^Ga! - X N (x = 0.33, 0.54) films were 
grown by MOVPE using the conventional method (the 
inset of Figure lc) to validate the prediction of the sur- 
face effect on Mg incorporation. As shown in Figure lc, 
the Mg concentration (C Mg ) on the surface is about one 
order higher than that of in the bulk for both samples. 
Although C Mg rapidly falls beneath the topmost surface 
(about 10 nm), C Mg is still several orders higher than the 
theoretical prediction by Equation 1. This phenomenon 
can be understood in terms of the competition between 
the Mg incorporation enhancement on the growing 
surface due to the surface effect and the Mg segregation as 
the epitaxy continues. Simply, when the surface with the 
enhanced Mg solubility is covered by newly added layers 
during further growth, most of these Mg segregates to the 
new surface to regain equilibrium because the surface 
transforms into a bulk. Meanwhile, considerable part of 
these Mg is frozen in because of solidification. As a result, 
the C Mg in the bulk is lower than that of in the final epi- 
layer surface but is much higher than the equilibrium 
value of the ideal bulk. Considering this competition, Mg 
incorporation can be modified by balancing the surface 
time and solidification time. 

As shown in Equation 2, the factor A/u also affects Mg 
incorporation. In principle, the chemical potential [i sen- 
sitively depends on the growth atmosphere (N-rich or 
Al/Ga-rich). The \a of a given species under equilibrium 
conditions is equal in all phases that are in contact [22]. 
Therefore, we can obtain 

f^Al/Ga = H-ANl/GaN ~ f*N 

= ^Al[bulk]/Ga[bulk] + f*N 2 + AH MN / GaN - ft N . (3) 

In addition, C Mg is limited by the formation of Mg 3 N 2 
to substitute Mg for Ga or Al as an acceptor [10]. This 
limitation meets the relation 



f*Mg^l/ 3 (f*Mg 3 N 2 - 2 f*N)- 



(4) 
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By substituting Equations 3 and 4 into Equation 2, we 
can obtain 



AH f >AE + ft AIN/GaN -l/3{i Mg3N2 -l/3{t 



(5) 



which, aside from A£, depends only on p Nf since the 
/^AiN/GaN and fi M g 3 N 2 are constants [25]. ^ N should be 
limited between ^ N (Al/ Ga-rich) < ^ < ^ N (N-rich) [11], 
namely, max(% 2 + AH AIN , p N2 + AH GaN )< (A N < (a N2 , to 
drive the source materials to form Al^Ga! _ alloys 
instead of the undesirable phases (bulk Ga, Al, and N 2 ). 
Our calculated AH GaN value of -1.01 eV is higher than 
the AH A i N value of -2.97 eV, which are consistent with 
the experimental values of -1.08 and -3.13 eV [25]. 
Therefore, as the growth condition varies from Ga-rich 
to N-rich conditions, changes from \a n<i + AHc aN 
to [t N2 . Thus, A/Ovaries over a range corresponding to 
l/3A// GaN of 0.337 eV, as shown in Figure 2a. This 
variation indicates that the N-rich growth atmosphere 
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Figure 2 Formation enthalpy difference of Mg Ga /Mg A | and C Mg 
profile of Alo.4 9 Ga 0 .5iN film, (a) Formation enthalpy difference of 
Mg Ga and Mg A , between Ga-rich and N-rich condition, (b) C Mg profile 
of Alo.4 9 Ga 0 .5iN film with three different Cp 2 Mg flows grown by the 
MSE technique. The inset in (b) illustrates the source supply sequence of 
the MSE technique, an ultimate V/III ratio condition is shortly produced 
during the interruption. 



favor the Mg incorporation effectively in AlGaN. Gener- 
ally, the N-rich condition is modulated by increasing the 
V/III ratio. However, for a fixed III flow, the Al^Ga! _ X N 
growth has an optimal V/III ratio for the best crystal qual- 
ity [13-16]. Nonetheless, the max flow limitation of the 
MOVPE system does not allow the V flow to be increased 
infinitely. Accounting for these limitations, an inspiration 
can be obtained from Figure lc, in which the protecting 
atmosphere with NH 3 flow just provides an ultimate V/III 
ratio condition (extremely N-rich) for C Mg enhancement 
when the epitaxy ends with the III flow becoming zero. 
Simultaneously, the stopped growth avoids the formation 
of low-quality Al^Ga! _ crystal. If this special condition 
is introduced as an intentional interruption during the 
continuous ^-Al^Ga! _ ^N growth, then the overall Mg 
incorporation could be improved. 

To validate this hypothesis, a growth interruption ex- 
periment was designed, as shown schematically in the 
inset of Figure 2b. We closed the metal flows (TMA1, 
TMGa, and Cp 2 Mg flows) three times. In these three 
periods (35 nm thick), different Cp 2 Mg flows (0.45, 
0.81, and 0.99 nmol/min) were applied to investigate 
the interruption effect systematically. Figure 2b shows 
the SIMS C Mg profile of Al 0 49 Ga 0 .5iN film across three 
periods. Obviously, C Mg exhibits three distinct peaks at 
the interruptions, and the profile of each interruption is 
similar to that of the final epilayer surface (Figure lc). 
Since there are always some Mg floating on the surface 
during growth because of segregation [26], the interrup- 
tion will drive the floating Mg to incorporate into the 
Al^Ga! _ ^N crystal, thus greatly enhancing Mg solubility. 
This result confirms that the Mg incorporation on the 
growing surface can be transiently enhanced further by an 
extremely N-rich condition interruption, thereby increas- 
ing the C Mg that would reside at the interrupting region. 

However, the C Mg enhancement at the interruption re- 
gion is much smaller than that on the final epilayer surface 
(Figure lc), and the C Mg far from the interruption region 
remains low. This result is caused by the wide interval be- 
tween consecutive interruptions, considerably decreasing 
the C Mg at the interruption regions and resulting in the 
non-uniformity of the C Mg distribution by Mg segregation 
and diffusion after interruption (Figure 3a). Therefore, the 
interruption interval, interruption time, and growth rate 
should play critical roles in affecting the C Mg overlap. As 
illustrated in Figure 3b, we further proposed the MSE 
technique, optimizing the interruption conditions, to in- 
corporate surface Mg atoms before they can re-segregate 
to the surface, thus further increasing the average Mg in- 
corporation and approaching a uniform Mg distribution 
over the entire AlGaN epilayer instead of being distributed 
locally. 

Three Mg-doped Al^Ga! _ X N (* = 0.54, 0.76, 0.99) 
samples were grown by using the MSE technique (the 
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Figure 3 Schematic diagram of the Mg incorporation behavior 
in the AIGaN grown by the MSE technique. As the interruption 
interval is long, only some peaks distribute locally at the 
interruptions after Mg segregation and diffusion (a), optimizing the 
interruption interval, a high and uniform Mg distribution over the 
entire AIGaN epilayer could be achieved (b). 



inset of Figure 2b). An optimized 2-nm interruption 
interval combining with 2-s interruption time were used 
for all samples, with Cp 2 Mg flux of 0.81 nmol/min. As 
shown in Figure 4a, the samples with different Al contents 
exhibit high C Mg range from 4 x 10 19 cm -3 to 5 x 10 19 cm -3 
and homogeneous distribution at a wide region as expected, 
whereas the C Mg of the samples grown via conventional 
method decrease with increasing Al content, which is 
consistent with the theoretical prediction. By compari- 
son, the average C Mg in the samples with different Al 
contents increase several times, and the enhancement 
ratios increase as the Al content increases, as shown in 
Figure 4b. Particularly, the enhancement ratio is approxi- 
mately up to 5 in the Al 0 99 Ga 0 .oiN. These results indicate 
that a high C Mg can be easily achieved in Al-rich AIGaN 
by combining the surface effect with the N-rich growth 
atmosphere modulation. 

High Mg doping was reported to result in Mg-rich 
precipitates. The primary Mg-rich precipitates were pre- 
sumed to be Mg 3 N 2 [27,28], which can be formed when 
Mg do not incorporate as acceptors in the desired substi- 
tutional sites. The substitutional Mg was suggested to be 
usually passivated by H during growth, and the corre- 
sponding Mg acceptor can be activated by postgrowth 
thermal annealing to dissociate the Mg - H complex [29] . 
The correlation between the substitutional Mg and H was 
verified by previous theoretical and experimental investi- 
gations [30,31]. Thus, the H concentration is most likely 
associated with C Mg if Mg is effectively incorporated in 
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the desired substitutional sites. The enhancement ratios of 
H concentration for the MSE technique increase from 1.2 
to 10 with increasing Al content, compared with that of 
the conventional method, as shown in Figure 4b. This 
simultaneous enhancement in H concentration demon- 
strates that the Mg was effectively incorporated in the 
desired substitutional sites by the MSE technique. In 
this work, the high C Mg is the important basis for im- 
proving the hole concentration in ^>-type AIGaN epi- 
layer. Besides the solubility limit, the high activation 
energy of Mg acceptors is another contribution for the 
low j?-type doping of Al^Ga! _ JM, leading to a low 



Zheng et al. Nanoscale Research Letters 2014, 9:40 
http://www.nanoscalereslett.eom/content/9/1/40 



Page 6 of 7 



acceptor activation probability [5,8]. In order to in- 
crease the overall p-type doping, more efforts on acti- 
vating the obtained high C Mg will be included in future 
progress. 

Conclusions 

The MSE technique, which utilizes periodical inter- 
ruptions under an extremely N-rich atmosphere, was 
proposed to enhance Mg incorporation, base on the first- 
principles total energy calculations. During the interrup- 
tion, metal flows were closed to produce an ultimate V/III 
ratio condition without affecting the AlGaN growth. By 
optimizing the interruption conditions, we obtained a high 
concentration and uniform distribution Mg in the AlGaN 
epilayer. The C Mg enhancements increase with increasing 
Al content through this method. Particularly, for the 
Alo.99Gao.01N, the enhancement ratio can be achieved up 
to about 5 and the final Mg concentration was determined 
to be 5 x 10 19 cm -3 . Meanwhile, the simultaneous increase 
of the H concentration confirms the Mg effective in- 
corporation in the desired substitutional sites instead of 
forming Mg 3 N 2 . The proposed approach, which is con- 
venient as well as effective, could be used as a general 
strategy to promote dopant incorporation in wide bandgap 
semiconductors with stringent dopant solubility limits. 
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